Aspects of Duality in Nodal Liquids 



Nick E. Mavromatos and Sarben Sarkar 
Department of Physics, Theoretical Physics 
King's College London 
Strand, London WC2R 2LS 

Starting from a microscopic t — J like model and a SU (2) spin-charge separation ansatz, a rela- 
tivistic continuum gauge lagrangian is obtained in the vicinity of a nodal point of the Fermi surface. 

■ The excitations in the pseudogap phase are described by topological excitations in the dual model 
' which has a Z2 global symmetry due to the effect of instantons. Confinement of spinon and holons 

emerge from this picture. The adjoint and fundamental strings are associated with stripes. As the 

■ spin gap decreases a local Z2 symmetry emerges. 
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I. INTRODUCTION 



High temperature superconductors (HTS) are materials which are believed by many to involve strong electron 
correlations [1] and so they are not in the universality class of BCS superconductors. Some believe that new concepts 
familiar in field theory such as strings arc relevant to their description. The phase diagram for HTS is complex 
involving antiferromagnetism, non-Fermi liquid behaviour, pseudogaps etc. [2] However, partly because the models 
are not solvable the basic mechanisms are still controversial. At the same time it is believed (or hoped) that the 
1 I gross features should have some underlying explanation based, not on the microscopic details, but on simple general 
principles which underpin the physics of all parts of the phase diagram. We shall present a model which has a plausible 
justification from one of the generic models that are in favour with the practitioners of high T c and is based on the 
inspired surmise that spin and charge somehow separate. This is familiar in condensed matter from studies of exactly 
solvable, typically lattice, models in one dimension. There is a sizeable and beautiful literature on this which comes 
under the heading of the Bethe ansatz, since in the distant past Bethe solved the one dimensional quantum Heisenberg 
model by ingenious methods that bear his name. With the passage of time the deep reason for Bethe's success was 
uncovered and this was linked to an infinite number of conservation laws in such models. The bold extension of these 
ideas to two dimensions has been made, reinforced by the presence of quasi-one-dimensional structures that will be 
' mentioned later. It is clear that the formal and rigorous grounds for Bethe's success do not hold in two dimensions; so 
we need to have a different route to the understanding of excitations in HTS. The spin-charge separation that we will 
t-H ' use does not directly spring from the Bethe ansatz result concerning the true excitations in one-dimensional systems. 
Rather it is a basis in which to define a model [3] . The originators had hoped that the basis would be a good starting 
point for the study of lattice models based on the experiences in one dimension. Irrespective of the pros and cons 
of such an approach we will employ it to derive a starting hamiltonian from a lattice t — J model. ( t denotes the 
magnitude of the hopping matrix element and J denotes the Heisenberg antiferromagnetic exchange.) Just as in the 
non-linear a model, there are constraints to make the dynamics more strongly coupled. The lattice constraint requires 
. at most one hole on a site. It has been argued that a simple one-band model (i.e. one involving one type of charged 
carrier 'orbital') may emerge as an approximation from a more realistic model involving Cu 3d x 2_ y orbital and the O 
(2p x , 2p y ) orbitals [4]. Even if the holes are localised on the O atoms, the lowest energy states can be (Zhang-Rice) 
O singlets formed with the hole on the Cu atom which can hop as a fermionic entity from plaquette to plaquette formed 
by nearest neighbour O atoms [5]. This is not a rigorous picture though, and features such as direct hopping of O 
holes, of magnitude t p , is ignored ( as opposed to that induced by hybridization of the O and Cu holes). The details 
are involved. Moreover the bare values are not necessarily the ones that are pertinent to the low energy physics (in the 
sense of running couplings used for the renormalisation group); we shall, as is traditional, ignore such complications 
and consider just the Zhang-Rice singlet picture. Given this input and the spin-charge separation ansatz, a certain 
amount of interesting and non-trivial phenomenology will arise. Inherent to the spin-charge ansatze are constraints 
and hidden gauge symmetries [3] [6] [7] . The exploitation of this gauge symmetry, using duality arguments developed 
for the description of confinement in 2 + 1 gauge symmetries, will form the central plank of our approach to stripe 
phases [8] and the resonant valence bond (RVB) states of high T c materials. In the dual description we shall see that, 
in general, a Z2 gauge theory will emerge (but from a perspective quite different to that used in some recent work 

Although various widely accepted effective lattice hamiltonians are simple to write down, they cannot be solved 
rigorously for planar systems. Consequently phenomenological effective continuum low energy hamiltonians, which 
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incorporate important symmetries, have been a focus of attention [10] . A spur to this has been recent experimental 
progress in establishing properties of the pseudogap phase in the underdoped HTS [11]. Whereas the superconducting 
phase has a well defined Fermi surface and quasiparticles, the underdoped region has a Fermi surface but not well 
defined quasiparticles. Moreover the pseudogap [12] itself has the <i-wave symmetry of the superconducting gap. 
There is some smoothness in the transition from the superconductor to the pseudogap. In the light of this, the low 
energy properties of the theory in the neighbourhood of the zeros or nodes of the superconducting gap may be a 
good starting point for investigating the physics of the pseudogap phase [10]. Near the nodes the theory can be 
described by a relativistic theory [7]. This is an example of the intimate relationship between different parts of the 
phase diagram. Moreover the inhomegeneities of the charge distribution inherent in the formation of stripes gives a 
physical separation of charges and spins. Our approach, which is restricted to low doping, will allow us to describe 
this separation in terms of the confining properties of 2 + 1 dimensional spontaneously broken SU (2) gauge theories 
[13] in a dual theory for nodal liquids. Moreover the existence of skyrmion textures centered on the doped holes [14] 
in this picture support an antiphase boundary condition for the staggered spins on either side of the stripe. 

Gauge theories have figured in theories of high temperature superconductivity since the early days of the subject [3]. 
Furthermore the strong coupling Hubbard model at half-filling, which leads to the Hciscnberg model for a quantum 
antifcrromagnet, was mapped onto a lattice gauge theory with SU (2) gauge group and Nambu-Dirac fermions [15]. 
Holes in antifcrromagnets were studied by Shankar [16] in one-space dimensions in the limit of large spin and small 
doping. This was generalised by Dorey and Mavromatos to two-space dimensions [6] . The large spin limit of the anti- 
ferromagnet is the same limit used by Haldane in his derivation of the a(CP 1 )-modcl [17]. Dorey and Mavromatos, 
and Shankar modelled the structure of copper oxide layers in HTS using a bipartite lattice. In the large spin limit 
holes hop only within each sublattice. This is due to assumed underlying Neel order and orthogonality properties 
of spin coherent states. The labelling of sublattices provides a flavour label for the fermions. The concepts of large 
spin and low doping will also be assumed in our analysis. However properties of the low energy theory around the 
nodes of the superconducting gap were not of interest at that time. It is interesting to note that, quite recently, by 
investigating Marshall wave functions for spin i Heisenberg antiferromagnets, qualitatively similar Berry phases were 
found for a hole performing a closed circuit [18]. Indeed these authors argued that this phase was fundamental to 
vanishing spectral weight and non-Fermi liquid behaviour. These phases can be incorporated in a statistical Us (1) 
gauge factor in the hopping factor of the Zhang- Rice singlet. 

In this article we will show how, starting from a microscopic lattice gauge theory, for which arguments have been 
given at length elsewhere [7] , some of the features of the pseudogap and stripes can emerge due to general features of 
symmetry and confining properties of gauge theories in 2 + 1 dimensions. Our aim is to incorporate certain symmetries 
within a spin-charge separation ansatz as a relativistic theory in the neighbourhood of a node. The framework is a 
development of the earlier work of Dorey and Mavromatos [6] and Farakos and Mavromatos [7] . The structure of the 
article is the following: in section II we review how one derives the continuum field theory of the nodal excitation 
under the non-Abelian spin charge separation ansatz in case of doped antifcrromagnctic systems with an approximate 
particle-hole symmetry. In section III we first derive in detail the dual theory in a Landau-Ginzburg approach in 
the case where the spinon degrees of freedom are assumed very massive so as to decouple. In section IV we discuss 
domain wall structures in this approach, which we connect with stripes in the original microscopic theory. In section 
V we incorporate the spinon degrees of freedom and show how, as the spin gap decreases, a local emerges in the 
theory. We also discuss antiphase properties of the staggered magnetization in the microscopic picture (as opposed 
to the dual one) by making the observation that the domain wall structure in the dual Lagrangian corresponds to a 
domain wall structure in the flux of the statistical photon, and stripe conductivity by discussing the Dirac equation 
of the nodal holons in such flux backgrounds. Conclusions and outlook are presented in section VI. 



It was noted by Affleck et al [15] that the spin operator at site i, in terms of the Pauli spin matrices a, can be 
written as 



where X = ( f t ) an< ^ Cl an< ^ ° 2 are ^ ne annihilation operators for up and down holes (which will actually 
\ c 2 ~ c i / 



represent the Zhang- Rice singlets). We will make a spin-charge separation ansatz [7] which has slave fermion character 
and incorporates the CP 1 constraint known from large-spin analysis at half-filling ( the Heisenberg limit). Our ansatz 



II. FROM MICROSCOPIC MODELS TO CONTINUUM RELATIVISTIC FIELD THEORY 
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X=(i\^)( Z l7 2 )=*Z (2) 



is 

r\ -r] J \ z 2 : 

where ip a and z a , a — 1,2 are fermions and bosons representing holes and spins repectively . The a index labels 
the sublattice structure of the bipartite lattice (relevant for the large spin analysis). The ansatz is valid in the low 
doping situation [7] and on assuming canonical (anti-)commutation relations for (ip a ) z a , canonical anticommutation 
for c a hold provided 

2 

ipi,iip2,i= = i>l ti ip\ ti , (^.^/M + = !• (3) 

13=1 

These are single occupancy constraints which implies that the ansatz requires strong coupling. There is also a 
certain redundancy in the description due to the following SU (2) gauge symmetry 

tf < -» Villi, Z, -» h\Zi where h, e SU (2) . (4) 
The global SU (2) spin manifests itself as 

Zi -» (5) 

where /i is a global transformation in SU (2). However, there is also a dynamical U s (l) gauge symmetry acting only 
on the V fields, which is due to phase frustration from holes moving in a spin background. Some arguments to justify 
this from a microscopic point of view have been given in reference [6] [18]. Consequently, this symmetry is associated 
with exotic statistics of the pertinent excitations [7] , which is an exclusive feature of the planar spatial geometry. 

By similar reasoning to Farakos and Mavromatos [7] we arrive at a generalised Hartree-Fock hamiltonian Hhf of 
the holon-spinon excitations: 

fa Ir { A 2 z}Ri jZj + h.c. J' (b) 

where Uij (Rij) is the U s (1) (SU (2)) link element between nearest neighbour sites i and j ( A\ and A 2 are constants 
resulting from assumed frozen amplitudes assumed typically in lattice gauge treatments). 

A necessary condition for a relativistic theory is that the fermions are coupled to the U s (1) gauge potential which 
fluctuates around a background with a flux of tt per plaquctte in each sublattice. The flux phase background [19], 
implying that a product of background U s (1) gauge potentials around each sublattice plaquette equals (—1), is 
essential for yielding relativistic fermions in the continuum limit according to the general analysis of Burden and 
Burkitt [20] . Such considerations lead to the form of our effective starting lagrangian. Instead of the 2x2 matrix 
structure of the fermion fields V it is convenient [7] to work in terms of two Nambu 2-component spinors: 

*I = (Vi,4) (7) 

and 

H = (*h,-4)- (8) 

In the background field of the flux phase of the U s (1) field a Kawamoto-Smit transformation [21] to another set of 

2-component spinor S Q is made. The components of S are linearly related to the components of V . In terms of these 
new spinors we obtain, in the action for the fermions, the following term which is necessary to have a continuum limit 
with Dirac structure: 

Sf = \ K ' Yl ( Sl ( _7 ^ Ri, fJ U ittl 'E i+fl + ^i+^^Rl jj}^ (9) 
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Here 5 denotes ET^o and the representation of the Dirac matrices that we use is 7^ = i (173, u\, cr 2 ) , /J, = 0,1,2. Ui.^ 
is the gauge link variable £/, i+ M . The gauge invariant kinetic terms for the link variables are taken to be 

S G = J2[p2(l-TrR p )+(3 1 (l-TrU p )} (10) 

p 

where p denotes plaquettes, the (3j are the inverse couplings; R p and U p symbolize the plaquette product over the 
SU (2) and U s (1) link variables repectively. The inverse coupling /3 2 is assumed large in order to have a conventional 
continuum Dirac form for the holon excitations. It should be noted that the kinetic terms for the link variables appear 
in any case if the spinon degrees of freedom are in the phase where they are massive (i.e pseudogap phase) due to 
quantum fluctuations . In terms of the microscopic parameters of the underlying t — J model, the SU (2) dimensionful 
coupling e (/3 2 oc has a spinon contribution of the form e 2 <~ J (1 — 5) where 5 denotes doping concentration [6]. 
Clearly the coupling due to the spinon contribution decreases as the doping increases. This implies that the gauge 
coupling for low doping is large. The kinetic term for the U s (1) gauge is not generated by quantum fluctuations of 
the spinons (which do not couple directly to the U s (1) gauge potential). This is equivalent to the condition /3\ ~ 0. 
Since /?i — 0, the Us {I) gauge field on the lattice can be integrated exactly [22], in the path integral for the partition 
function, to leave 
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VRVEVE exp(—S e ff) (11) 
where 

Seff =f32^2( 1 - TrR p) + J2 l °S h (VV^) (12) 

p i,n 

and 

y itl = -k 2 Tt (m« (- 7 „) R^M^^A,) ■ (13) 

k is related to k' and is a constant depending on the microscopic energy parameters of the underlying theory and 
in particular the doping. M denotes a composite variable M a b. a f3 — S& i( gS 0)a , which can be expanded in terms of a 
complete basis of bilinear fields [7] : 

AfW = ]T $ a (i) a a + S 4 (i) 1 2 + i ((S,) 4 (i) r + ( F ^a (0 l^a j • (14) 

a=l \ a=l ) 

Here 

$1 = -i (SiS 2 - S 2 Si) , $2 = (HiS 2 + S 2 Si) , $ 3 = (S1S1 - S 2 S 2 ) (15) 

form an adjoint representation of SU (2), 

(A ll ) 1 = i (Si<7 M S 2 - 220-^21) , (A fl ) 2 = (Sia M S 2 + S 2 cr p Si) , (-4. M ) 3 = (Sict^Si - S 2 ct aj S 2 ) (16) 

form a vector adjoint representation, the S4 and («S M ) are two singlets and ctq = —ias,a\ = a\,oi = cr 2 . Since 5 is 
grassmanian, the Taylor series in S truncates and so 

- log /* r (2^/y-) - - w „ + iy 2 , - + (17) 

which is an exact expression. This effective potential allows us to consider the dynamical formation of various 
condensates of fermionic bilinears. It is also necessary to pass from the path integral over fermionic variables S, H 
to the bosonic variables in terms of which M is written . The transformation has a non-trivial Jacobian [22], which 
leads to additional terms in the effective action of the form — ^ ^ logdet MjM i+A1 (where i is a lattice site index 
and fi is a space-time index ). In the ground state a scalar condensate along the 0-3 direction is formed, i.e. 

(M^) = u<r 3 + (18) 

which implies that [7] [22] 

($ 3 ) = u. (19) 



4 



This state is a mimimum of the effective potential. On writing 

Ri„ = cos (| B it , |) + ia.B ¥ sin (j B l{1 |) / | | (20) 

we find that 

XogPj (2^/y-) ~ Ml ({B^y + (B^) 2 ) + interaction terms (21) 

with M| = k 2 u 2 . Moreover a parity-invariant mass is also generated for the fermions by the strongly-coupled Us (1) 
interactions in G, entailing a gap for charge excitations (consistent with the pseudogap phase). The SU (2) group is 
broken down to U (1) and there remains only one massless vector boson, B^, the dual photon. An important aspect 
of this U (1) symmetry is that it is compact and allows for monopole solutions which are interpreted as instantons 
in 2 + 1 dimensions. This will lead to a complete breaking of the symmetry of the dual theory and also a small 
non-perturbative mass for the dual photon consistent with the pseudogap phase. On taking a naive continuum limit, 
a low-energy gauge theory results, which will be the starting point of our analysis. For a nodal liquid there should 
be a Lorentz invariant ground state and so it is to be expected that only the scalar multiplet $ a , a = 1, - • • ,3 
can have non-zero expectation values. Lorentz invariance requires that the ground state expectation values of F a ,n 
should vanish. In addition the Vafa-Witten theorem [23] rules against the formation of parity violating states (i.e. 
non-zero values of condensates of Si). For the properties of the nodal liquid we shall assume, for these reasons, that 
^,54 and do not appear in the effective theory. They will be important sufficiently far from the nodes. The 
fields that will appear in our continuum theory will be A all ,& a and z a in such a way that SU (2) gauge symmetry 
is manifest. Higher order couplings, from power counting, are expected to be less dominant for small momenta, and 
consequently should not be important for an effective low energy theory. The U s (1) field will not explicitly appear 
but its fundamental effect of spontaneous breaking will be incorporated through a Higgs potential term. Similarly 
the fundamental constraint in the theory, because our formulation is for low doping, can be subsumed into another 
Higgs potential term with suitably chosen coefficients to implement the constraint as tightly as required. Coalescing 
all these considerations allows us to define our starting continuum effective lagrangian [24]: 

C M = --F^F^ + \ (D^V) 2 MV 2 - A (^ 2 ) 2 + \ (d^z^ + £ {? if (22) 

where D^ ab = d^5 a b + ee cab A c ^, D^ ab = d^S a b — ^eA^a^ b . This theory involves fields in both adjoint and 
fundamental representations of SU (2). The former is due to the charge degrees of freedom and the latter is due 
to the spin degrees of freedom. This will be important for our future considerations. The above lagrangian captures 
the qualitative features of our approach, i.e. the symmetries and the mechanism for their breaking. For [i 2 < the 
intuitive classical analysis indicates spontaneous symmetry breaking, and, for ri 2 > 0, full symmetry. A rigorous 
analysis indicates no breakdown of analyticity in going from positive to negative /j 2 . Because of symmetry-breaking, 
compact U (1) gauge theory emerges out of an explicitly SU (2) symmetric theory, and this is also necessary in the 
gauge formulation of the CP 1 model which requires a U (1) gauge field [25]. This emergence of a compact U (1) is an 
important difference from the earlier work of Dorey and Mavromatos [6] . 

III. DUALITY 

We will first consider the z fields to be very massive (i.e. a large spin gap) and so for the time-being ignore them. 
The spin gap disappears where superconductivity starts up [26] and so, we would expect that as we approach the 
superconducting transition, our arguments will need to be modified. The dynamical symmetry breaking crucial to our 
analysis is induced by the holon condensate $ fields. By working in the unitary gauge (953) 7^ 0, the massless gauge 
field is A^ which we will call the statistical photon (as opposed to the electromagnetic photon). The remainder of 
the gauge fields have a mass rmjv and the Higgs field has a mass m#. These conclusions are based on a perturbative 
analysis. However there is also a crucial non-perturbative effect which is due to instantons [27] which gives the 
statistical photon a mass m^-The theory in the absence of the z fields is the well-known Georgi-Glashow model which 
has been a toy model for the study of confinement based on Zn vortices. However in our case it has quite a different 
interpretation and, for the group SU (2), is not a toy model. We will identify some of its non-perturbative features in 
terms of the pseudogap phase of HTS. The addition of the z degrees of freedom changes a discrete global symmetry to 
a local one [28]. Since we are considering zero temperature, fields are defined over 3-dimensional space. Consequently 
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the monopoles of the cuclidean 3-dimcnsional version of the theory can be interpreted as instantons in one lower 
spatial dimension. These monopoles will affect the disorder fields used in discussions of duality. At its root duality 
is an elaboration [29] of the one due to Kramers and Wannier [30] for the Ising model. A phcnomenological study 
of nodal liquids using duality, linking the XY model and the Landau- Ginzburg model for a scalar field, has been 
pursued by Balcnts, Fisher and Nayak [10] but is in quite a different spirit. 



We will assemble a dual description bit by bit. If we ignore the effects of monopoles, absent in a non-compact U (1) 
theory, then our starting theory, represented by a lagrangian Cah , is the Abelian-Higgs {AH) model. The associated 
hamiltonian density H, in the temporal gauge, takes the form: 



H = \ {B 2 + El) + ttV + | D^\ 2 - M 2 H 2 + A | 



where [i and A are constants, Di4> — {di + ieAi 



i = 1, 2 and ir — dnd, w' = do 



Aj (t, x^ , Ei {t, ~af ) = 5ij5 (if — x/^j and (t, x^ , n {t, ~a?) = 5 (it — . 
solutions which are given asymptotically as 



(23) 



the non-trivial commutators are 
This model has classical vortex 



i arg("x*) 



a?) — ► (p e 

|"af | >oo 



Ai (t, it) — > — d t arg {it 

\t?\ — >oo e 



(24) 
(25) 



where arg {x) is the angle between x and one of the co-ordinate axes. The quantum version of the vortex is represented 
through the vortex creation operator [31] [32], the disorder variable V, which is dual to Ai [33]. Its commutation 
relations with A; and <t> are 



V {t, If) A, {t, t) = (Ai (t, t) - -d\ arg ( f - T?)J V {t, a?) 



(26) 



and 



V- 1 {y) 4> {x) V {y) = e iai **-*4> (f) . (27) 
The equal-time commutation relations (26) and (27) can be satisfied with the following representation of V 

{Xj - Vj) 



V {If) = exp I - / d z y 



{ x - y ) 



Ei {Y) + arg ("a? - 1?) j { it) 



(28) 



In fact, in the absence of fields in the fundamental representation, i.e. the z fields in the lagrangian, V {Hf) is local 
[24], [34]. As the z fields become lighter and so need to be taken into consideration, V will become non-local and the 
description has to change in a profound way. From (26) it follows that 



2tt 

[V (t, a?) , B {t, It)] = ~—V {t, H?)5{x-y) = -gV {t, !?)8{x- y) 



(29) 



and 



[V{t, !?),<S>{t)} = -gV{t, If) (30) 

where $ (i) = J d 2 yB (t, 17*) , B is the magnetic field and g is the 'magnetic' charge conjugate to e. Hence V can 
be seen as creating a fundamental unit g of magnetic flux which is compatible with being a disorder variable. This 
will be a key observation for us. For the Higg's phase the flux symmetry is not broken and (V) = whereas in the 
Coulomb phase (V) ^ [34]. The photon is the Goldstone boson of the broken symmetry. The phase transition can 
be interpreted as a condensation of vortices. We should not get confused here by the mention of phase transitions 
where previously we stated that there is no phase transition. There is no contradiction since we are supressing the z 
degrees of freedom. It is the presence of the fundamental fields which implies the lack of a phase transition. 
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A Landau-Ginzburg approach to the low-energy theory, initiated by 't Hooft [31], leads to a dual lagrangian d 
which has the form 



C d = d lt V& t V-Sl(V*V) (31) 

where f2 is a Higgs-like potential. The field V is not restricted to have strictly unit magnitude since this allows some 
of the features of quantum fluctuations to be seen at the classical level. This is an effective dual picture valid near the 
cross-over from Higgs to confining behaviour ( since a phase transition is absent) . The classical solitons in this dual 
description are the elementary excitations in the original quantum theory [31] [24], [35]. Before we consider the form 
of the solitons, it is necessary to incorporate the effect of instantons (i.e. monopoles). The magnetic current is no 
longer conserved but satisfies 

dpF" (x) = \d^ v °F va (x) = k (x) (32) 
where k (x) is the non-zero magnetic charge density ( in analogy to the U (1) charge current ). ft is given by 



k (x) = —£ abc £ ^d^ (x) d v 0> (x) d a (p c (x) 



(33) 



and [36] 



= F^ a - -e abc tp a D^ b D^ 



(34) 



where (p a 



--. Consequently, when there exist monopole solutions, the magnetic flux is no longer conserved. 



This complication results in a form for V which is identical to that already given with Ei (1/*) replaced by F° l (I/*). 
For a monopole centred at the origin of 3-dimcnsional cuclidcan space 



~ Air 

= — <5 3 (x) . 



(35) 



On integrating over t ,t Q x R 2 and remembering that there is a mass gap we conclude from (35) that 



$ (to) -*(t ) (36) 

This explicitly shows that the flux is not conserved. In the presence of a finite number of monopoles and anti- 
monopoles, the flux can change with time by a multiple of ^2L. The unitary operator U a (t) for a flux transformation 
is given in general by 

U a (t) = e lea ^ (37) 
but a is restricted. In fact for instanton configurations with Kronecker index n, the operators 

U n = e** n * (38) 
labelled by an integer n is an invariance of the Hilbert space provided 

Un (to) U- 1 (to) = e 2 ™ lfc = 1 (39) 

where k is an integer. Moreover the commutator (30) still holds (with the appropriate forms of the operators) [24]. 
From (30) the values of the flux on physical states can only be multiples of ^ . Hence only the operators U n with 
n = and n = 1 are independent. Consequently because of instantons the U (1) flux symmetry reduces to a Z2 
symmetry and so Cd has to have non-[7 (1) symmetric terms which nontheless are Z 2 symmetric. Hence £4 becomes 
[31] 

C d = dpVdpV - A (V*V - £ 2 ) 2 - ^ (V* 2 + V 2 ) + C Ova d u V*d x V) 2 . (40) 
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The four derivative term in Cd is important for describing the breaking of the adjoint string which is a higher energy 
phenomena and so is dominant at short distances. The system has a mass gap since the phase of V develops a mass 
2\fh. This is interpreted as the photon acquiring a mass because of instanton effects and is connected with the absence 
of long range phase coherence in the pseudogap phase [7] . However there is no genuine phase transition between the 
confining and Higgs 'phases', and so Cd which is most easily derived through semi-classical methods in the Higgs 
phase (weak coupling) is valid also in the confining phase. The difference lies in the spectral properties i.e. whether 
rriA <C ma ( in the Higgs phase) or vice versa (in the confining phase). The mass tha is associated with the phase 
of V and the mass mn with the radial part of V. From semi-classical and weak coupling calculations the parameters 
in Cd can be identified in terms of the masses and couplings of the microscopic lagrangian Cm as follows: 



t2 



e 



2 



~ _ 2n 2 m 2 H 

* ~ 2 

h = rriA 

C0C ^- (41) 



IV. ADJOINT TEXTURE 



There are various topological textures which are possible in the dual theory [24], [35]. Quantum features in the 
original theory can be seen at a classical level in terms of these textures. The winding number of the V field in the 
texture is given by the charge of the U (1) symmetry in Cm- This is encoded in the expression for the U (1) current: 

J, = -i^e^ p d v (V**dV) . (42) 

Clearly the charge is non- vanishing only where the phase of V is singular. 

An isolated vortex defect has an energy which diverges linearly with the size of the system. The hedgehog configu- 
ration is quadratically divergent and is not energetically favourable owing to the breakdown of the Z7" (1) symmetry 
to Z 2 (due to monopole effects in the Higgs phase). In the confining phase the same symmetry picture arises due to 
infrared singularities [37] which is to be expected if there is no genuine phase transition. This defect is the adjoint 
string [35] and corresponds to the charged SU (2) gauge particles. A finite string and also finite energy would arise 
from a Zi vortex and anti- vortex. The particles that these dual structures represent are bound states of the statistical 
gauge field. A statistical gauge field cannot be seen directly but it will leave its imprint on the properties of the holons 
and spinons. We will describe this in the remainder of this article. Its thickness can be estimated to be inversely 
proportional to itia [24]. In the confinement phase this is small and the string is narrow. Let us recall that e 2 ~ J 
and so we are in the strong coupling situation for low doping. We contend that the adjoint string is the basis of stripe 
behaviour. Firstly let us discuss how long this might be. Typically stripes are long and thin. We can write V as pe l6 . 
For the string to end it needs a vortex structure of V at which to terminate. The size of a vortex core is the region for 
which p differs from its vacuum value. Moreover at short distances the dominant energy term is £ (e fj. v \d v V* d\V) 2 . 

Now the size of the core is O ^^-^and so the scale of the core energy £ can be estimated by (mjjp 4 . For weak 

coupling we can get more information since we can identify the dual couplings with those from the microscopic theory 
( cf (41)). Although we are really interested in strong coupling, the absence of phase transitions emboldens us to use 
weak coupling and extrapolate. We will not be accurate but the hope is that the qualitative inferences will hold for 
strong coupling. Hence we have that the core energy density £ satisfies £ ~ (m^e 4 . It is known that the adjoint 
string tension a a ~ e 2 rriA [24] and so the length of the stripe L before it snaps, i.e. the typical length of a stripe, is 
given by 

o A L ~ {m%e 4 ( — ) . (43) 

This gives L <~ -For thin stripes it is necessary that L » It may not be possible to establish without a 

detailed numerical analysis. However given the absence of a true phase transition (in the full theory) we might be 
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able to glean something from an extrapolation of the weak coupling analysis. From the semi-classical analysis of 
Polyakov [27] 



^ 16^r< ^ T W £ / m*\l 

where e is a slowly varying function but limited in range to 1 < e < 1.787. This exponential factor when substitute 
into the estimate for L gives an indication that for strong coupling the ratio of the length to the width can be a large 
factor. Owing to the uncertainties of the extrapolation we cannot be more quantitiative. In the strong coupling limit 
with wla > run the width of the stripe is inversely proportional to tua- In the absence of the pseudogap the width 
of the stripe goes to infinity since to a — * and so the stripe structure disappears. These one dimensional structures 
which are related to the adjoint field lead to one-dimensional charge structures. These will play a catalytic role for 
the formation of domain structures. There are however other crucial features which need to be incorporated within 
this picture. The first is that holon transport will be along the stripe. The second is the antiphase nature of the spin 
structure where the staggered magnetisation will point in opposite directions on either side of the stripe. In order to 
investigate these aspects it will be necessary to consider the spinon degrees of freedom which are in the fundamental 
representation. 



V. SCALAR FUNDAMENTAL 



The model of (22) has additional terms in the scalar field z which transforms as the fundamental representation of 
SU (2). This actually is necessary for establishing that there is an analytic path in coupling space which connects the 
Higgs and Coulomb 'phases' [38]. We have relied on this result in the previous section even though the z-degrees of 
freedom did not appear directly.The dynamical symmetry breaking of SU (2) to U (1) allows the CP\ model for the 
z s which is necessary for a description of the magnetic degrees of freedom. The L of 22 has a global U (1) symmetry 
given by 

and the associated quantum number will be called a 'magnetic' number since this symmetry is just part of our global 
rotation symmetry (5). In our effective dual lagrangian this global U (1) will play an important role. 

The previously introduced vortex field V is neutral under this transformation and is associated with the degree of 
freedom derived from the frustration of hole motion in an antifcrromagnetic background. In the prescence of the z's 
V is non-local [28]. This is related to the fact that the Wilson loop has a perimeter law independent of confinment or 
the | Higgs phase and so V could not be a local order parameter. Once V is allowed to be a local order parameter 
it is possible to give a standard and heuristic argument by splitting the minimal surface inside the loop into little 
areas; an area or perimeter law that differentiates between confinement and the Higgs phase ensues. We will find 
that there are new solitons for V , the fundamental string, in the presence of these magnetic degrees of freedom. 
At energy scales where magnetic numbers appear, in the dual representation it is necessary to add a magnetically 
charged complex field U to incorporate this phase symmetry. In our effective lagrangian the field will be important for 
inducing spontaneous symmetry breaking of the U (1) gauge group down to Z 2 - Moreover defects (vortices) of the U 
field ( given by the zeros of U) will be contained within the core of the V defect and represent spinons. The U field is 
thus a dual represenation of the z field. By duality vortices in z ( i.e. skyrmions) will represent smooth configurations 
of U ( or U condensates). Since in the pseudogap phase there is a spin gap, it is necessary for the spinon fields to be 
massive (with mass M z ). In our formalism this is seen from the constraint on the z fields. In the presence of matter 
fields in the fundamental representation, the hamiltonian involving V fields is not globally Z 2 invariant but has a 
local Z 2 invariance. Consequently it is necessary to have a lagrangian which is manifestly Z 2 gauge invariant. The 
amalgamation of the above considerations leads to a dual lagrangian Ld [28] with the following structure: 

+ l\(d»-ib»)U\ 2 -\ v (V*V — fj 2 ) 2 

-\ u (U*U - u 2 ) 2 - £ (V 2 U* + V* 2 U) + ( {e^ x d v V*d x V) 2 (46) 

where f^ v = d^b v — d v b^ . We consider X u >• X v as well asu>/j. The gauge transformations of U and V are 
interlocked because, when V — ► Ve la , it is necessary that U — » Ue 2ia (manifest from the term proportional to £). This 
lagrangian is a generalization of that given earlier (40). Although we cannot provide a rigorous justification for the 
effective lagrangian Ld , we will show that it embodies the important symmetries that are present in Cm (22). 
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It will first be shown how in the the absence of fundamental fields Ld reduces to the lagrangian in (40). We first 
note that the following reparametrisation [28] of e,u and z, with the assumption that y, x and k are independent of 

leads to the desired reduction. In the limit M z — > oo, e — > and the vector boson mass mj, which satisfies m 2 = e 2 u 2 , 
remains finite. This results in a lagrangian with terms which coincide with 40 together with terms involving decoupled 
W fields. The above parametrisation can be deduced for large M z with y oc A 2 as will be shown below. 



The lagrangian Ld has two global U (1) symmetries. The associated currents are 



e 

and 



j" = (48) 



/„=<W- (49) 
In terms of Cm we can make the following identification: 

^7, - if (50) 

where the current iff connected with the transformation 45 is 

if = i(z a '^^z a ). (51) 
Also the charge Q = J j°d 2 x can be expressed as 



= J d 2 xdi ^4 



Q = / d 2 xdi -reafj (52) 



and so is related to the vorticity of the magnetic current. We will now demonstrate that this vorticity of the 
magnetic current is created by the operator V. Using Gauss's law in 2-dimensions, we can rewrite the operator (28) 
as the exponential of a line integral over a path C x which is the branch cut of the multi-valued function arg (y — x) at 

x. For clarity we will denote this form of the operator by Vc x (= cxp J c dsiSijEj^j). This operator is identified 

with V in Ld after gauge fixing. It can then be shown that [28] 

V^jt 1 (a) V Cx = if (*) + nnf (x) 6 (x e C) z a * (x) z a (x) (53) 
where nf (x) is the normal to C x at x. The vorticity density w of the magnetic current is defined as 

- = (f^f-) ■ ( 54 ) 

From arguments similar to that used in (53) 

V^u>(v)Vc.=w{y)+Tr6<V(x-y) (55) 

and so Vc x creates vorticity it, i.e. half a unit of vorticity. From the gauge couplings of U and V in L it is clear 
that U will create twice the vorticity produced by V. 

This derivation allows us to identify e 2 in Ld- Indeed, from the relations (48) and (49) we can deduce that 

- = (J ft ) ^ 
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and so e 2 ocz" z b . In the effective theory it is a good approximation to replace z z by its expectation value and 

so e 2 oc (z b * z b ) . 

This identification allows us to verify the scaling of e in (47). Indeed, after performing the frequency integration, 
the leading estimate for (z b z b ) is 



f 2 1 A 2 
z - z" ) = / eTp oc — (57) 

which allows us to infer that y cx A 2 where A is a suitable high momentum cut-off that defines the effective theory 
and so verify (47). 

We now need to examine the influence of the massive spinon fields. For energies much less than the scale u 2 we 
saw that the U degree of freedom was frozen and that Cd (40) was recovered which led to the adjoint string. However 
there are other configurations which we have so far neglected and are crucial to the development of the domain 
wall structure central to our considerations. For these there are points where U vanishes and so the unitary gauge 
transforming U to u cannot be fixed. The vortex of U carries magnetic quantum number. From Lr>, by examining 

the covariant derivative structure, it is readily seen that for finite energy for both V and U vortices it is necessary 
that a single winding of the V field is accompanied by a double winding of the U field. Moreover, since sufficiently 
far from the U vortex a unitary gauge form of U is valid, the field V forms a one-dimensional adjoint string . In 
the absence of U vortices it must be accompanied by an anti-defect which, for Z 2 , is the same as the defect itself. 
This is the incipient charge order which, on introduction of the z 's will lead to a special type of spin ordering of the 
staggered magnetization. Furthermore charge transport will also be along the texture. Before that it is necessary to 
establish the domain wall structure in V in the presence of fundamental charges, through an argument due to Fosco 
and Kovner [28]. We recall that the fundamental charges in the dual picture are represented as the vortices in U. 

Let us see how this happens by using the parametrisation U = ue l8 and V — pe l{p . The dual lagrangian Cd has 
the form 

C D = l -^pf+ p 2 (d^ -h^j j - 4Cp 2 {e^xVpfrip) 2 - X v (p 2 - p 2 ) 2 - 2£,p 2 u cos (2<p - 9) (58) 
and at long wavelengths the equation of motion for V reduces to 

d 2 p = p (d^ip - - AX V (p 2 -p 2 )p- ^pu cos (2ip - 6) (59) 



and 



[p 2 [d^~\h] ) =4£p 2 usin(2^-- 0). (00) 



From (58) it is clear that winding numbers in U and are in correspondence while V winds only half as much, i.e. 
far from a vortex ip = |. This is the crucial intuitive reason for the formation of the domain wall for V . Consider an 
adjoint string with vortex and antivortex at its ends and then introduce a spinon at each end. We would like investigate 
the effect this has on the nature of the texture for V. Explicitly consider U vortices, the spinons, in the plane at 
x = (±a, 0). Far from the vortices 

U(x) = uexp(i6(x)) (61) 

with 

6 (x) = tan" 1 ( -^-) + tan" 1 ( -^-) (62) 
\xi-aj \xi+aj 

due to the winding of the two vortices. As x 2 — > oo, 9 (x) — ► 7r and (p — > -|. Similarly as x 2 — > — oo, 9 (x) — ► — it 
and ip — > — j. Consequently there is a domain wall formed. 



11 



We will now discuss the antiphase properties of the staggered magnetiztion in the microscopic picture ( as opposed 
to the dual picture). The above domain wall structure translates into a domain wall structure in the flux of the 
statistical field due to (30). We need to find solutions of the CP 1 field z in the presence of this domain structure. 
The configuration taken up by this field is that of vortices on one side of the domain wall and an anti- vortices on the 
other. We can infer this from the form of a single static CP\ vortex 

,^ / COS^fi e -W*) \ 

sin ii e l^) ) (63) 

which is centred at the origin and where r = \x | , x = ^,arg(x) = arctan (§7) an d f (r) = 2 arctan^, s being 
an arbitrary scale factor which will be the thickness of the domain wall i.e. s ~ ^j. ( In terms of the staggered 

magnetisation n (n, = z a * {<Ji) ab z b ) this configuration is that of a skyrmion viz. n{x) = (sin/(r) x, cos/ (r)).) The 
winding number of the vortex is the flux of the 'magnetic' field B. The form of B is 

s sin ( 2 arctan 



b & = act? ^/Trv ™ 

y/x( +xi, (xj +x% + s 2 ) 

and so falls off rapidly away from the vortex. Consequently on one side of the domain wall there is a skyrmion and 
on the other an anti-skyrmion in terms of the staggered magnetisation. Because of the lack of rotational symmetry 
in the presence of the domain structure, we expect, in keeping with the symmetry of the domain wall (or stripe), a 
row of vortices and anti-vortices on cither side of the stripe. This means that there is an anti-phase structure for the 
staggered magnetisation. 

In the microscopic picture the nodal fermions satisfy a Dirac equation (with real time) which near the domain wall 
takes the form 

(i Jd - e a e /A - e a m) ip a = (65) 

where ei = 1 and e 2 = — 1- (Because of the signs associated with the masses there is no parity breaking in the model 
[6].) A related construction in a non-relativistic theory for staggered magnetic fields can be found in [39]. We consider 
a domain wall along the x-axis with A (x) = B (0, \x 2 \ , 0). On making the ansatz for ip of a propagating solution 

ip a {xo,Xi,x 2 ) = exp (i (luxq ~ kxi)) Xa (^2) , (64) 

where Xa is a 2-component spinor for each a we deduce that 



|^-f +0- 1 -e a sgn(x 2 )^ Xq = 



(65) 



where 9 = u 2e B an( ^ £° = y ~eE + e a e B \x 2 \) from (65). Although these equations hold near the stripe, we can 
check whether there are propagating solutions localised in x 2 . If so then the stripes could conduct. In our situation 
£1 — > 00 and £2 — > — 00 when \x 2 \ — > 00. For s=lwc impose the boundary conditions at £1 ^00 and y = 0; for a = 2 
we impose the boundary conditions at £2 — * —00 and y = 0. The symmetric (+) and antisymmetric (— ) solutions can 
be written as 

« (s > - ( ) » > »<»> - ± ( d d :% ) » < «. 

in terms of parabolic cylinder functions Dg and similarly the solutions vanishing at £ 2 = — 00 are 

My) - ( ° e{ 7lt ) v>o, My) = ± ( \} { J^ ) y < o. (67) 



We are still left with one boundary condition at y — 0: for each sublattice a, we impose the condition 



lim Xa(y) = lim Xa(y), (68) 

y~*0+ y^0- 
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such that we have for both sublattices 



De-i(e a tio) = ±De(e a t ) , £ = k^j A (69) 

which gives us the dispersion relation 0(£o) or equivalently u)(k). It is possible to solve these relations which are 
consistent (and in the non-relativistic limit give rise to modified Landau levels). The leading asymptotic behaviour 

Dp (z) ~ e~~z p , |^| — > oo and |argz| < -7r (70) 

and so the wavefunctions are localised on the scale of the 'magnetic' length and so localised conduction along the 
stripe is consistent within our framework. 



VI. DISCUSSIONS 



In this work we have started from a SU(2) <g> U s (l) gauge field theory of spin-charge separation, based on an 
approximate particle-hole symmetric formulation upon doping lightly from half-filling. The strongly coupled U s (l) 
group is a symmetry only of the fermionic sector, expressing appropriate frustrations of holons. Once U s (l) is 
integrated out in the path integral one arrives at an effective lagrangian with a broken SU(2) — > 17(1) phase. In this 
picture, the presence of doping is responsible for a dynamical breaking of the initial SU(2) spin symmetry. 

Non-peturbative effects, due to the compactness of the U(l) unbroken subgroup, are responsible for giving a mass 
in the associated gauge boson ("statistical photon"). This leads to a pseudogap phase, that is a holon mass-gap phase 
without phase coherence. In this article we have discussed in some detail some properties of a formalism dual to 
the above mentioned effective lagrangian, which proved very useful in shedding light on various important physical 
properties of the original theory. In particular, we have argued in detail how the dual formalism can explain in a 
natural way the stripe phase of the underlying microscopic theory. We have associated such stripes with appropriate 
topological textures (domain walls) of the dual lagrangian. The role of vortices of the dual lagrangian as appropriate 
dual configurations of the original spinon (magnon z) degrees of freedom has also been pointed out. Moreover, anti- 
phase properties of the staggered magnetization in the microscopic picture (rather than the dual one) have been 
studied by making the observation that the domain wall structure in the dual lagrangian corresponds to a domain 
wall structure in the flux of the statistical photon. 

We would now like to conclude by making a few speculative remarks on the role of nodal excitations to the passage 
from the pseudogap to the superconducting phase of the underlying antiferromagnetic system. The important point 
to realize is that, in the context of our relativistic nodal theory, superconducitivity arises in the way explained in [6] 
[7], only in the case of a strictly massless statistical photon, i.e. a photon that remains massless non-perturbatively. 
This photon plays the role of the massless pole in the current-current correlator, which is the basic Landau criterion 
for superconductivity. In this model, the pseudogap phase -studied in the present paper- can be succeeded by a 
superconducting phase if and only if the non-perturbative monopole effects, responsible for the statistical photon 
mass, are suppressed. 

We have conjectured in [7,22] that this may happen in the case where there is a dynamical super symmetry in the 
effective theory between spinon and holon degrees of freedom. Such a situation has been argued by Mavromatos and 
Sarkar [40] to characterise specific points in phase space of some extended t — j models under the spin-charge non- 
Abelian separation ansatz (2). Subsequently, it was argued by Alexandre et al [41] that the continuum composite field 
theory of such supersymmetric theories, obtained after the integration of the strongly coupled U s (l) group, exhibits 
extended N = 2 supersymmetry, due to the low-dimensionality; this arguably can provide the necessary mechanism 
for the masslessness of the statistical photon, and thus a passage from the pseudogap to the superconducting phase. 

It should be remarked that the presence of extended supersymmetries opens up the way for some exact results 
in the phase diagram of such systems. An interesting question is to formulate the duals of such supersymmetric 
theories and study their properties in detail along the lines presented above, e.g. the fate of textures etc. This is 
left for future work. It remains to be seen, of course, whether realistic microscopic models for antiferromagnets can 
be constructed which, in some regions of their parameter space, could exhibit dynamical supersymmetries between 
spinon and holon degrees of freedom, capable of explaining phenomenologically the rich phase diagrams of the high- 
temperature superconductors observed experimentally. We believe that the present work, along with those in [40] 
[41], constitute useful contributions to this programme. 
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